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^ : ABSTRACT 

^ ■ We investigate the properties of highly compressible turbulence, the compressibility 

■ arising from a small effective polytropic exponent 7e due to cooling. In the limit of 
O ! small 7e, the density jump at shocks is shown to be of the order of , much larger 
^ \ than the jump associated with high-Mach number flows in the isothermal regime. In 

the absence of self-gravity, the density structures arising in the moderately compressible 
• case consist mostly of patches separated by shocks and behaving like waves, while in 

Q ! the highly compressible case clearly defined long-lived object-like clouds emerge. The 

^ \ transition from wave-like to object-like behavior requires a change in the relative phase 

^ ■ of the density and velocity fields analogous to that in the development of an instabil- 

ity. When the forcing in the momentum equation is purely compressible, the rotational 
energy decays monotonically in time, indicating that the vortex-stretching term is not 
^ \ efficient in transferring energy to rotational modes. This property may be at the ori- 

■ gin of the low amount of rotation found in interstellar clouds. Vorticity production is 
found to rely heavily on the presence of additional terms in the equations, such as the 
Coriolis force at large scales and the Lorentz force at small scales in the interstellar 
medium, or on the presence of local sources of heating. In the presence of self-gravity, 
we suggest that turbulence can produce bound structures for 7e < 2(1 — 1/n), where 
n is the typical dimensionality of the turbulent compressions. We support this result 
by means of numerical simulations in which, for sufficiently small 7e, small-scale turbu- 
lent density fluctuations eventually collapse even though the medium is globally stable. 
This result is preserved in the presence of a magnetic field for supercritical mass-to-flux 
ratios. At larger polytropic exponents, turbulence alone is not capable of producing 
bound structures, and collapse can only occur when the medium is globally unstable. 
This mechanism is a plausible candidate for the differentiation between primordial and 
present-day stellar-cluster formation and for the low efficiency of star formation. Finally, 
we discuss models of the interstellar medium at the kpc scale including rotation, which 
restores a high-7e behavior. 
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Introduction 



Vstrophysical flows are expected to be highly rur- 



bulent, with a wide variety of mostly compressible 
energy sources at various scales. Thus we can expect 
the nature of astrophysical turbulence to be quite dif- 
ferent from standard laboratory turbulence, which is 
mainly a vortical phenomenon. If the flow is sub- 
ject to strong heating processes and radiative cooling, 
with thermal time scales substantially shorter than 
the dynamical ones ( Spitzer fc Savedoff 1950 ), as is 
the case in the interstellar medium (ISM) or the post- 
recombination medium in cosmology, rapid thermal 
equilibrium can be established, and the fluid behaves 
as a polytropic gas in which P oc p^", where 7c is an 
equivalent polytropic exponent. Such systems have 



been discussed both in the linear (Elmegreen 1991) 
and nonlinear regimes (Vazquez-Semadeni et al. 1995, 
hereafter referred to as Paper I; Passot et al. 1995, 
hereafter Paper II). In particular, in the models of 



the ISM at the 10-1000 pc scales we have presented 
in previous work, the thermal time scales range from 
10"^ to 10~^ times the dynamical timescales (Paper I; 



see also Elmegreen 1993a). For standard cooling hmc- 



tio ins (e.g., Dalgarno fc McCray 1972 ), which depend 



he mctallicity and ionization fraction among other 
factors, 7e can take values smaller than one and even 
negative (leading to the isobaric mode of the thermal 
instability). As 7c approaches zero, pressure gradients 
are 



almost absent, leading to a rather inelastic behav- 



ior With a very high compressibility, a behavior which 
is reminiscent of Burgers flows even when the Mach 
number based on the isothermal sound speed (the 
quantity normally reported in observational studies) 
is of order unity. 

In this paper we investigate the influence of the 
degree of compressibility of the flow (measured by ei- 
ther the effective polytropic index 7c or the actual 
Mach number Ma) on both the properties of the tur- 
bulence and the mechanisms of formation of density 
structures. It is important to predict whether tur- 
bulence is dominated by vortical or compressible mo- 
tions, depending on the type of forcing imposed on 
the flow, and the additional physical processes con- 
sidered. In particular, it is of interest to investi- 
gate the nature of the transfer between compress- 



ibla and rotational modes of the velocity field. Ihis 



tional to compressible modes, and only in an efficient 
manner above Ma ~ 0.3 (see also Passot & Pou 



quet 1987). In the present paper we extend their re- 



sults to the highly compressible regime, and use two- 
and three-dimensional numerical simulations to show 
that, for purely compressible forcing (such as that 
originating from stellar winds), vorticity decreases 
steadily. We also discuss alternative vorticity pro- 
duction mechanisms, such as the baroclinic term or 
the presence of shear, the Coriolis force, and mag- 
netic fields. Throughout the paper, we address the 
discussion mostly towards processes in the interstel- 
lar medium at various scales, although we consider 
generalizations which should be applicable to a wider 
variety of astrophysical flows, in particular in Cosmol- 
ogy and in molecular clouds. 

In the context of the ISM, both the amplitude and 
persistence of density fluctuations arc highly depen- 
dent on the degree of compressibility of the flow, and 
thus the latter can affect the cloud and star forma- 
tion mechanisms. In the linear regime, Elmegreen 
(1989) has proposed that the gravitationally unstable 
scales are greatly reduced in the presence of dissipa- 
tion through cloud collisions, while Renard & Chieze 



(1993) have studied the modiflcations to the Jeans 
length in the vicinity of the thermal instability. In 
addition, various workers have discussed the forma- 
tion of self-gravitating clouds induced by external ve- 
locity fields (e.g 



Hunter & Fleck 1982; Hunter et 



al. 1986| ; [Elmegreen 19934 Paper | |Padoan 1995| ) 



problem has been investigated m the low-to- moderate 
cornpressibility regimes by Kida & Orszag (lyyua,b), 
who found that transfer occurs mostly from rota- 



although without consideration of the direct depen- 
dence on an effective polytropic index 7e. In this 
paper we also discuss, based on numerical simula- 
tions, the ability of turbulence to form gravitation- 
ally bound structures which are the sites of star for- 
mation, as well as their characteristic scales, in re- 
lation to 7c. This problem will be addressed in the 
context of both a polytropic gas (either with a sin- 
gle value of 7c or a piecewise index depending on the 
density) and a fully thermodynamic model of the ISM 
including parameterized heating and cooling, as in- 
troduced in Paper I and, with rotation and magnetic 
field ( Paper II ). Note that in this paper we only con- 
sider the large-scale effects of turbulence in the forma- 
tion of density structures, neglecting its well known 
effects for cloud support at small scale s (e.g., Chan-| 
drasekhar 1951; Bonazzola et al. 1987; Leorat et al 



199C; Vazquez-Semadeni & Gazol 1995; see also the 
review by Scalo 1987| and references therein). How- 



ever, our results should also be applicable to a cloud- 
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fluid in which the pressure includes a sub-grid tur- 
bulent component (e.g., iHunter & Fleck 1982|; IScalo 



& ptruck-MarceU 1984| ; |Elmegreen 199lD ; sub-grid 
cloud collisions provide an effective power-law "cool- 
ing" mechanism (e.g., Elmegreen 1991). 

The outline of the paper is as follows. In the 
next section we describe the numerical approach and 
define the various models. Section || is devoted to 
the influence of various degrees of compressibility on 
randomly forced compressible turbulence, when the 
forcing is only in the velocity equation and is only 
potential, and to the production of vorticity in the 
highly compressible regime by the vortex-stretching 
term and other mechanisms, in two and three dimen- 
sions. The next Section then deals with the impli- 
cations for scale-free astrophysical flows when self- 
gravity and magnetic fields are included. Section ^ 
considers the results of the previous sections in the 
context of models of the ISM at t he kpc scale based 
on the one presented in Paper II . Finally, § || is the 
conclusion. 

2. The models 

In this paper we use several sets of hydrodynamic 
equations to simulate highly compressible flows, rang- 
ing from a simple barotropic flow to the full model of 
the ISM used in Paper II (with the addition of winds 
or a random forcing Fc in the velocity equation). The 
latter reduces to the following equations, which rep- 
resent the ISM on the plane of the galactic disk: 



dp 

di 



^ + V.(pu) 



(1) 



- (V X B) X B - 2f7 X u - t/8 V^u 
P 

V2{V^u+\w • u) + Fs + r 



de V^e 

— -I- u • Ve = —(7 — l)eV • u + kt 

ot p 



(2) 



dt 



Fd + r, - pA, 

V X (u X B) - i/gV^^B + r/V^B, 



1, 



(3) 

(4) 
(5) 



P = (7 - l)pe, (6) 
rd(x,<)=ro(p/pic)-", (7) 

{constant if p(x, to) > Per 
and < t - to < Ms , (8) 
otherwise 

and 

A = A,r^- for T, <r<r,+i, (9) 



where 














Ti 


^ 100 


Ai 


= 1.14 X 10^5 


Pi 


= 2 




T2 


= 2000 


A2 


5.08 X 10^6 


P2 


= 1 


5 


T3 


= 8000 


A3 


= 2.35 X 10" 


P3 


= 2 


867 


Ti 


= 10^ 


A4 


= 9.03 X 10^8 


Pi 




0.65 




= 4 X 10^. 













The temperature T is measured in K and A in erg s~^ 
g~^ cm^. Note that the stellar heating Fg is spread 
out on a few pixels around a stellar center by a con- 
volution with a Gaussian profile. In eq. (||), is the 
time at which the density threshold pcr is reached at 
point x (see Paper I and II for further details). 

As usual, p is the density, u is the fluid velocity, 
e is the internal energy per unit mass, P is the ther- 
mal pressure, B is the magnetic induction, O is the 
angular velocity of the rotation, and cj) is the grav- 
itational potential. The ratio 7 = Cp/c^ of specific 
heats at constant pressure and volume respectively is 
7 = 5/3. The temperature is related to the internal 
energy by e = c^T . The nondimensionalization is the 
same as in Paper II; the resulting nondimensional pa- 
rameters are J, the number of Jeans lengths in the 
integration box, whose side is 27r, and Ma, the Mach 
number corresponding to the characteristic velocity 
Mo at the initial mean temperature of the flow (note 
that we take Cv = [7(7 — l)Ma]~^)- The magnetic 
field is written as B = BqGx + b where Bq represents 
the uniform azimutal component of the field, and b 
a superimposed fluctuating component. The shearing 
forcing Fg is accomplished by fixing one component 
of the largest Fourier mode in order to impose a ve- 
locity profile Ux = —0.375 siny; the additional term 
Fc is discussed further below. 

These equations are solved both in two and three 
dimensions with periodic boundary conditions, using 
a Fourier pseudospectral technique with a resolution 
of 128^ collocation points for two-dimensional runs 
and 64"^ for the three-dimensional ones. A hyper- 
viscosity scheme with a V® operator is used, which 
confines viscous effects to the smallest resolved scales 
and allows for the use of much smaller values of the 
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kinematic viscosity V2 and magnetic difFusivity r\ than 
what would be possible when i/g = 0. The adjunc- 
tion of the second-order viscosity and difFusivity, com- 



pared to the model used in Paper I and Paper I], 



has been made necessary in the highly compressible 
regime in order to "filter" out the oscillations in the 



vicinity of strong shocks (Passot & Pouquct 1988) 



In addition, a mass diffusion /iV^p is added to the 
continuity equation in order to smooth out the den- 
sity gradients, thus allowing the simulations to reach 
higher r.m.s. Mach numbers. The values of the vari- 
ous coefficients for the 128^ runs are = 8 x 10""'^^, 
V2 = ri = 2x 10"■^ ^ = 3 X 10^2 and Ma = 1 through- 
out the paper. For the 3D runs we take — 10~^, 
V2 = 2.5 X 10~^, and ^ = 4 x 10"^. We refer the 
reader to Paper I] for further details concerning the 
technique and the model terms. The additional ac- 
celeration Fc introduced in the velocity equation has 
been taken either as a random vector field with a k~'^ 
power-law spectrum peaked at a wavenumber fcfor, 
white noise in time and with an adjustable ratio of 
compressible to solenoidal components, or as a wind- 
mimicking momentum source of the form Fc — — V?/;, 
where il> is proportional to the stellar heating Pg. Fi- 
nally, the initial conditions for all variables are Gaus- 
sian fluctuations with random phases. For the ve- 
locity field, the fraction of solenoidal to compressible 
energy is an additional adjustable parameter. 



According to the results of Paper I and Paper I], 
due to the short thermal time scales, the flow be- 
haves essentially as a piecewise poly tropic gas, except 
in the regions of star formation. As a consequence, we 
also solve a simplifled set of equations (hereafter de- 
noted ppm, for Piecewise Polytropic Model) in which 
we eliminate the internal energy equation and use a 
piecewise polytropic equation of state given by 



P 



ap 



,7c i 



P n9' 



A; 



7^2 



(10) 



with i — 1^ 5, in conjunction with the five temper- 
ature domains of eq. (^. Note that in the barotropic 
case the nondimensionalization of the equations leads 
to an explicit factor of Af^^ in the pressure, the sound 
speed in the equilibrium state of temperature Toq be- 
ing given by Cg = (7e/7M^)Teq. The effective poly- 
tropic index is given, in terms of the exponents Pi in 
(||) and a in eq. (|^), by: 



7e 



1 



(1 + a)M 



(11) 



present notation stresses its dependence on the inter- 
val of temperature in the cooling function. Finally, 
the boundaries of the different temperature ranges in 
the cooling function for the full model correspond to 
densities given by 



P 1 
' Ai ' 



Pi 

rp l + o 



(12) 



Note that 7o,j was noted 7e// in Paper II, but the 



For convenience, we adjust the degree of compress- 
ibility of the flow (the value of 7eJ by varying the 
exponent a of the diffuse heating. This is equiva- 
lent to varying the exponents j3i in all intervals of the 
cooling function. Although this is not a very realis- 
tic approach, it provides us with the simplest means 
for varying the flow's "hardness", i.e. the flow's re- 
sistance to compression. Actual modifications to the 
cooling function will be discussed in a future paper 
( [Franco ct al. 1996| ). 

A further level of simplification, useful to single 
out purely fluid dynamical properties of supersonic 
turbulence, consists of using a single polytropic expo- 
nent 7c. In this case, the pressure is taken as P = 
(7cM^)~^p'''=, leading to a sound speed Cc = (I/Mq) 
in the equilibrium state p = 1. In the following, we re- 
fer to these equations as the B-model, for Barotropic 
model. 

Throughout the paper, we will refer to a rather 
large number of simulations, whose most relevant pa- 
rameters are summarized in Table 1. 

3. Some properties of highly compressible 
flows 



3.1. Density structures 

As discussed in the introduction, low values of 7o 
lead to very high compressibilities, and in fact the 
degree of compression reached in this case is much 
higher than the one obtained by increasing the Mach 
number. This can be seen by calculating the density 
jump X = P2I P\ across a shock in a barotropic gas of 
index 7c. It is easy to show that X satisfies 

- (1 + 7cM2)A + 7cM2 ^ (j^3) 

where M is the Mach number upstream of the shock. 
From this equation, we recover the fact that the com- 
pression ratio for an isothermal shock (7c = 1) is M^, 
but we also see that, for 7c <C 1, A" is approximately 
given by X (1 -I- 7cAf2)^/^=. -p^j. ^ fixed value of 
M, this expression tends to e*^ as 7c — > 0, a density 
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jump which can be much larger than the isothermal 
one. 

We now discuss the nature of the density conden- 
sations present in the flow as observed in our simula- 
tions, and their relation to the value of the polytropic 
exponent. For simplicity, we consider the B-model 
and the ppm equations, with a purely compressible 
random forcing in the velocity equation. 

Even though the boundary between moderate and 
high compressibility is not sharp, the two regimes lead 
to quite different density structures. In the former 
case, the density field consists of patches separated by 
shocks, and whose lifetime is the sound crossing time 
within the patches. This is illustrated in fig. |l| (left), 
where we display a typical snapshot of the density 
field in a simulation labeled run 63, with 70 — 1.5 
and fcfor = 4. Note that the highest peaks form at the 
intersection of shocks. Such a behavior is analogous 
to what is observed in decaying flows, both in 2D 



dPassot et al. 1988[ ) and in 3D ( [Porter et al. 19941 ). 

As 7e is lowered, the density peaks become higher 
and narrower, and have longer lifetimes since the ef- 
fective sound speed decreases. This is illustrated in 
fig. |l| (right), showing the density field of run 65, with 
7o = 0.3. In fact, when using the ppm equations with 
the fiducial values of the parameters used in Paper II 
[a = 0.5), 7e can locally reach zero for intermediate 
densities, and we then observe a transition of behav- 
ior from wave-like to object-like. We exemplify this 
behavior in fig. ||, which shows the density field of 
two two-dimensional simulations using the ppm equa- 
tions, the first one with a = —0.5 (run 54, left), and 
the second with a = 0.5 (run 53, right). 

In order to understand the distinction between 
wave-like and object-like behavior, consider a passive 
scalar, which is defined as a quantity s obeying the 
pure advection equation ds/dt = — u ■ Vs, and con- 
sider its motion relative to a density perturbation. 
We define an object as a density perturbation within 
which the passive scalar remains for a substantial pe- 
riod of time, whereas in a wave the passive scalar and 
the density fluctuations drift apart. For example, a 
one-dimensional wave packet p — p{x — uot), moves as 
a whole at speed Uq. Adams et al. (1994) have shown 
that, in the context of a logatropic equation of state 
of the form P ~ c^p + Po log(p/ pn) such wave packets 
can behave as solitary waves. Conservation of mass 
leads to (u — uo)p = est, showing that the fluctua- 
ti ons of p a nd u are "in phase" (fig. ||; see also fig. 2 
of Paper Il| ) . The scalar s, advected at the flow veloc- 



ity u ^ Uq, will as a consequence drift with respect to 
the density perturbation. On the other hand, in the 
case of an object (i.e., a density clump), the velocity 
field is either convergent towards its center ("out of 
phase" with the density perturbation), or is zero in- 
side and contains two accretion shocks at its edges so 
that the passive scalar remains within the boundaries 
of the density fiuctuation. The transition from wave 
to object can spontaneously happen when the flow 
becomes unstable, e.g. thermally or gravitationally. 
In the presence of turbulence this behavior is contin- 
uously approximated as the polytropic exponent ap- 
proaches zero. In this case, the flow behaves almost 
as a Burgers flow, with large density peaks at shock 
locations instead of the density jumps appearing in 
flows with finite pressure. Therefore, the clouds have 
long lifetimes and are only affected by collisions with 
other clouds. Note that in the limiting case of zero 
pressure gradient, cloud (shock) collisions become en- 
tirely coalescing ("sticky"), since disruption requires 
the presence of vorticity, whose dynamics is discussed 
in the next section. 

3.2. Evolution of vorticity in highly compress- 
ible turbulence 

3.2.1. The role of vortex stretching 

It is first useful to write the equations for the po- 
tential vorticity Wp = uj/p = (V x u)//9 and the di- 
latation V • u in the barotropic case. Neglecting the 
dissipation terms proportional to i/g and V2 , they read 



dt 



u ■ Vwp = Wp • Vu 



(14) 



9V-U 



dt 



-u-VxtiJ- 



070 



7c 



- 1 



. 

Note that these equations are asymmetrical: dilata- 
tion can easily be produced out of rotational motions 



(from the 



-UJ 



term), whereas vorticity production 



can only originate from the vortex-stretchi ng term 



Vu in three d imensions (3D) (see also Kornrc- 



ich fc Scalo 1996 . In the two-dimensional inviscid 
case, it is well known that there is conservation of any 



functional of the potential vorticity (e.g., Pedlovsky 



1987) in the absence of shocks. When oblique shocks 



are present, the vorticity contained in them will be 
dissipated much more rapidly than elsewhere in the 
flow in the presence of viscosity. Moreover, since cjp 
obeys a passive scalar equation (in 2D), the vorticity 
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w (which is larger in dense regions) will tend to be 
confined inside clouds or at their edges in the highly 
compressible regime, in which, as mentioned in § [3.1| , 
clouds behave as objects. This was observed in the 



simulations of the full model with a — 0.5 (Paper II 


) , sumed 


but 


not in those of 


Paper I 


with a — 0, i.e. with a 1980|), 



higher 70. 

We now examine whether the vortex-stretching 
term can contribute significantly to the transfer of 
energy from the compressible to the solenoidal modes 
of the velocity. Note that, even though in 3D the pres- 
ence of the vortex stretching term does not allow us 
to draw any conclusion regarding the conservation of 
Up, it is nevertheless relevant to consider the equation 
for this quantity, since it avoids the need to take into 
account the compression term — cjV-u in the vorticity 
equation, which only reflects mass conservation. 
Xida fc Orszag (1990a) have investigated the trans- 



fer between internal modes and compressible and ro- 



tational kinetic modes in three dimensional, forced, 
compressible turbulence at moderate Mach numbers. 
They find a negligible transfer between the compress- 
ible and rotational kinetic modes. Here we extend this 
analysis to highly compressible flows, using the ppm 
equations. Starting with 100% solenoidal motions, 
and stirring the flow with a purely potential accelera- 
tion, we observe, for a wide range of values of 7e, that 
the ratio of rotational to total kinetic energy per unit 
mass [| es/eic always decays, as shown in fig. ^. To- 
gether with other runs discussed below in §3.3.1, this 
figure shows the evolution of eg/ek for two 3D runs 
at resolution 64^ solving the ppm equations, one with 
a = 0.5 (run 60), and the other with a — —0.5 (run 
67). The latter is thus "harder". Note also that, after 
the initial transients have subsided, the rate of decay 
appears to be independent of 7c, at least for the range 
of values considered here. Thus, the stretching term 
does not appear to contribute to energy transfer from 
the compressible to the rotational modes. 

The above result is rather unexpected, since vor- 
tex stretching is the main mechanism for vorticity 
production in incompressible turbulence. Instead, for 
highly compressible flows, vortex stretching appears 
to contribute mostly to the destruction of potential 
vorticity. We suggest the following mechanism. Let 
us first consider the case in which potential vorticity 
is a small perturbation on an irrotational flow. Since 



eq. (|lj) is identical to that for the magnetic field, it 
is legitimate, in analogy with the kinematic dynamo 
mechanism, to investigate whether cUp can increase on 
scales large compared to that of the potential flow, as- 



Radler 



namo happens through the presence of helicity in the 
flow. In the present case, however, the helicity of the 
basic flow is identically zero, rendering the growth of 
potential vorticity unlikely at dominant order. Never- 
theless, we must consider the possible ampliflcation of 
ujp on a scale comparable to that of the basic flow (as 
in the context of the fast dynamo problem) . In order 
to investigate this problem, a knowledge of the struc- 
ture of the basic flow is necessary. As discussed in the 
previous section (see fig. Q), the highly compressible 
regime is dominated by the presence of locally two- 
dimensional^^ shocks (see also Kadomtsev fc 
Petviashvili 1973|). As already observed by Kida & 



Orszag (1990b), vortex lines tend to be antiparallel 
with density (and velocity) gradients. A possible ex- 
planation of this observation is that, on both sides 
of a shock, the velocity gradient matrix has positive 
eigenvalues, which tend to amplify vortex fluctuations 
in the direction perpendicular to the shock, while the 
vortex stretching term would have no effect on vor- 
tex lines orthogonal to the density gradient. However, 
while the vorticity is amplifled in these expansion re- 
gions, it is also advected towards the shock from both 
sides, where it is then destroyed since the velocity 
gradient matrix has negative eigenvalues in these re- 
gions of strong compression. In summary, in highly 
compressible regimes, shocks are expected to act like 
sinks in which potential vorticity is drained from the 
flow. 

In the fully nonlinear case where rotational modes 
are initially of the same order of magnitude as com- 
pressive ones, the numerical experiments shown in flg. 
^indicate that the decay of rotational energy persists, 
suggesting that the above mechanism is still at play. 

3.2.2. The baroclinic term 

Another important source of vorticity in the fully 
thermodynamic case is the so-called baroclinic term, 
VP X Vp/p^, which is known to be active behind 



curved shocks or at th e ir collisions ( Passot fc Pou- 



^Note that the integrated energies per unit volume and per unit 
mass typically differ by no more than a few percent. 



quct 1987| ; [Fleck 199l| ; [Klein fc McKee 19941) . To 
test its effect, we have performed a new simulation, 
run 68, identical to run 60 described above, but us- 
ing the full thermodynamic equations. We flnd that 
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the rotational modes still decrease steadily in time, 
as also shown in fig. ^. This indicates that the baro- 
clinic term has a globally negligible contribution to 
the vorticity production, being overwhelmed by the 
draining effect of the vortex-stretching term in nega- 
tive dilatation regions. This remark is corroborated 
by the observation made by Porter et al. (1995)] that 
on average the baroclinic term is ~ 20 times smaller 
than the stretching term in compressible flows forced 
by an external shear. Note that, behind the intersec- 
tion of s hocks, vorticity creation can be locally im- 
portant ( Passot fc Pouquct 1987 ), whereas in other 
regions the press ure and density gradien ts are always 
close to parallel ( Kida fc Orszag 1990b| ). 

In fact, vorticity production behind shocks does 
not require the contribution of the baroclinic term 
( Hayes 1957] ), the production mechanism relying on 
sin iple geometrical considerations. In the present sim- 



ula tions such vorticity production is not seen, possi- 
bly due to the lack of resolution or to the fact that 
the shocks that form remain straight (occurence of 
such vorticity production has been observed in high- 
resolution simulations of collisions of strong SN shock 



waves with spherical molecular clouds (Klein & Mc 



Ke(! 1994)). The presence of the baroclinic term does 
not change this conclusion and is only important, as 
we shall see, when thermal heating is present (§ ^.3.3 ). 



3.3. Other mechanisms of vorticity produc- 
tion 

Since the nonlinear and baroclinic terms have been 
found in the previous sections to be globally insuffi- 
cient for vorticity production and maintenance in the 
presence of high compressibility and purely compress- 
ible forcing, it is important to examine alternative 
sources of vorticity in astrophysical flows. As an ex- 
ample, we will consider the case of the ISM at the kpc 
scale, within which possible agents for vorticity pro- 
duction are the Galactic disk rotation, various heat 
sources and the magnetic field. 

3.3.1. Coriolis force 

It is well known that the Coriolis force can con- 
vert converging motions into shearing ones. In fig. ^ 
we also show the time evolution of the ratio eg/ek for 
the three dimensional run labeled run 69, analogous 
to run 60, but with the addition of the Coriolis force 
with r2 = 0.4 in code units, cor respondin g to the fidu- 
cial value 27r/2 x 10^ yr~^ (see Paper II ) . It is notice- 



able that compressible and rotational modes almost 
reach equipartition after 4 x 10^ yr. Compared to the 
run without the Coriolis force (run 60), with the same 
random forcing, run 69 develops 30% less density fluc- 
tuations but twice as much kinetic energy, while the 
internal energy remains approximately at the same 
level. The increase of the level of kinetic energy is 
probably due to the fact that compressible motions 
dissipate faster than rotational ones and thus, for a 
given energy input rate, the equilibrium kinetic en- 
ergy will be larger if rotational modes are excited. 
Additionally, it is interesting to estimate the impor- 
tance of the Coriolis force at different scales. For econ- 
omy we do this in two-dimensions. The Rossby num- 
ber (measuring the relative size of the Coriolis force 
with respect to the nonlinear advection term), scales 
as L^/'^ if we assume that the velocity dispersion u 



scale s roughly as L^^"^ ( Larson 1981 ; Falgarone et al 
1992| ). A run using — 0.05 (run 85), a value ap- 



propriate for scales ^ 25 pc, ends up with only ~ 5% 
of rotational energy, while a run with = 0.4 (run 
57, appropriate for the kpc scale) typically contains 
~ 50% of rotational energy. Interestingly, this value 
is very similar to the 3D case, run 69. This indicates 
that the Galactic Coriolis force is likely to be an im- 
portant source of vortical motions at the kpc scale, 
but not at the molecular cloud scale. 

3.3.2. Large-scale shear from differential rotation 

Another agent due to the Galactic rotation is the 
shear arising from its differential character. A run 
with an imposed shear as performed in Paper 11 cor- 



responding roughly to Oort's A constant (run 59), 
shows a selective creation of vorticity of the same 
sign as that of the imposed shear. For this run, 
Gs/ek ~ 30% on average. 

3.3.3. Local heating 

The stellar heating Fg provides another source of 
vorticity in the flow. Indeed, at locations where ther- 
mal energy is injected, the baroclinic term will be 
important, since the gradients of density and tem- 
perature cease to be correlated. This effect is most 
noticeable in the dynamics of expanding shells, which 
propagate without deformation when only wind-type 
forcing is included in the equations, while they dis- 
tort and fragment in the presence of stellar heating. 
This is illustrated in fig. |[ which shows the density 
fields for runs 28 (left) and 35 (right). These runs 
start with purely compressible velocity fluctuations 
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and differ only in that the forcing in run 35 is in the 
velocity (through V-0; cf. § while run 28 has ther- 
mal energy forcing only (through Fa). The production 
of solenoidal energy is substantial in run 28 in which, 
during the period of shell breakup, eg/ek ~ 0.5, while 
in run 35, e^jey barely reaches 0.01. Note that the 
low value of Cg/ek in run 35 is in agreement with the 
result of j |3.2.2 that the production of vorticity behind 
curved shocks is negligible in the absence of external 
heating. 

Note that the breakup of shells in the presence of 
heating can also be interpreted as a consequence of 
lin ear Vishniac-typc instabilities, occurring in fronts 



bet ween thermal and ram pressures ( Vishniac 1983 \ 



However, the present simulations do not allow us to 
decide whether the shells break due to vorticity pro- 
duction by the baroclinic term, to the linear Visniac 
instability, or both. This is because in these runs ac- 
tually a thin shell of star-forming sites forms imme- 
diately behind the gas shell, causing the baroclinic 
term to be large there. One indication is that in 
run 35, with wind-type forcing, the nonlinear Vish- 



niac i nstability (Vishniac 1994; see also Hunter et al 
19^ ^tevens, Blondin fc Pollack 1992| ) does not de- 
velop, suggesting that the present simulations lack 
the necessary resolution. A detailed examination of 
this problem would require monitoring of the baro- 
clinic term and possibly higher-resolution simulations. 
This, however, is beyond the scope of the present pa- 
per. Here we just record the importance of local heat- 
ing as a source of vorticity. 

3.34. Magnetic field 

Finally, we discuss the effect of the Lorentz force 
on the transformation of compressive motions into ro- 
tational ones. For this purpose, we have performed 
four runs with the 'p'pm model labeled runs 86, 84, 87 
and 88, identical to run 53 except for the inclusion 
of a uniform initial magnetic field along the x-axis, 
with strengths of, respectively, 0.05, 0.3, 1, and 3, in 
nondimensional code units (at — \^ the Alfven 
velocity Va = uq, see Paper II). A small fluctuating 
field with an rms amplitude of 0.01 is also included. 
In fig. II we show the time evolution of the fraction of 
specific kinetic energy in solenoidal motions Cg/ek for 
the four runs. The fraction of energy in solenoidal 
modes is seen to increase with until equiparti- 
tion with compressible modes is reached at values of 
Bo > 1- The fluctuating magnetic energy o beys the 
same trend (fig. 0), as already mentioned in Paper II 



in the context of the full model. In fig. ^ temporal 
oscillations are observed, whose frequency is roughly 
proportional to B^. However, these oscillations do 
not appear to be due only to Alfven waves. Instead, 
fully developed MHD turbulence is present, as can be 
inferred from the presence of locally strong currents 
- indicating small-scale intermittent magnetic struc- 
tures - and a power-law magnetic spectrum (fig. H). 
The increase of the solenoidal energy in the form of 
fully developed turbulence with the uniform compo- 
nent of the magnetic field is remarkable and might, 
at first sig ht, appear contra dictory to common as- 
sump tions ( Elmegreen 1993b , sect. VII; Mouschovias 
199l| ). This phenomenon of creation of vorticity by 



magnetic fields is important as well in the context of 
accretion disks ( Vishniac fc Diamond 1992 ). 

The relative importance of the Lorentz force in the 
momentum equation increases as smaller scales are 
considered. For instance, taking typical values of the 
fiuctuating magnetic field strength of ~ 5/iG at the 
kpc scale and 30/iG at scales ~ 5 pc, and typical 
velocity dispersions of respectively 10 and 2 km s~^, 
the ratio of the Lorentz force to the nonlinear advec- 
tion term is roughly 20 times larger at the 5 pc scale. 
Thus, the production of vorticity from compressive 
motions due to the Lorentz force is most efficient at 
molecular cloud scales. 

4. Scale-free flows v^rith self-gravity 

4.1. Influence of 7c on the formation of col- 
lapsing structures 

The results presented in the previous sections have 
implications for star formation, since the ability of 
turbulence to form self-gravitating clouds depends on 
the compressibility of the flow. A simple estimate 
for the required values of 7e for this to occur can be 
given as follows. Consider a medium with an effective 
polytropic exponent 70, for which the effective Jeans 
length is given by 



in 



Gpo 



2-7c 



1/2 
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(16) 



where C\ is the isothermal sound speed such that 
P = cfp and Lj is the Jeans length based on Ci. The 
critical density at which the scale L becomes unstable 
is then given by pj oc l2/(7c-2)_ order to account 
for turbulent compressions acting on n directions, we 
consider a volume V = L"-Lq^", where L is the side 
of the volume which varies upon compression, and Lq 
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is the side that remains unahered. The critical mass 



various workers (e.g., Hunter & Fleck 1982; Bonaz- 



to destabilize this volume is thus given by 



zola et al. 1987; Lcorat ct al. 1990; Padoan 1995) and 



Mj oc L 



(17) 



If Afj is a decreasing function of i, turbulent com- 
pression can produce gravitationally unstable struc- 
tures, similarly to Hoyle's (1953) notion of hierarchi- 
cal gravitational fragmentation. Setting the exponent 
of L to zero then gives the critical value of the poly- 
tropic exponent as a function of the dimensionality of 
the compression n. In particular, 



7c 



for n 
for n 
for n 



(18) 



This result thus recovers the well-known critical val- 
ues of 7e for the existence of equilibrium solutions 
of self-gravitating polytropic gas configurations (e.g., 
Chandrasekhar 1961), although within the different 
context of n-dimensional turbulent compressions. In 
the present case, this criterion describes the threshold 
value 7c of 7e at which the internal energy of a turbu- 
lent fluctuation increases more slowly upon compres- 
sion than its gravitational energy. Conversely, above 
7c, a fluctuation can never be rendered unstable upon 
compression. Note that for 70 < 0, no compression 
is necessary because the medium is then thermally 
unstable. 

If we extend this result to non-integer values of 
n, accounting for clo uds with fractal boun daries of 
typical dimension 1.4 ( Falgaronc et al. 1991 ), and re- 
calling that the compression i s likely to b e driven by 
quasi-planar shocks (see also Fleck 1996 ), we obtain 
that < 7c < 1, which is probably not unrealistic. 

In order to examine the interplay between turbu- 
lence, compressibility and self-gravity, we have per- 
formed a series of forced runs with the B-model, us- 
ing barotropic exponents of 0.9, 0.3 and 0.1 (runs 106, 
107 and 108 respectively) and the same Jeans length 
ieff = 27r/ J (which is independent of 70, a fact follow- 
ing from the form of the pressure for the B-model), 
with J = 0.9. Note that all these runs are stable 
according to the linear criterion. The evolution of 
the maximum density within the integration box for 
runs 106, 107, and 108 is shown in fig. ^ The turbu- 
lent density fluctuations never grow in the case of the 
high-gamma run, but in the case of runs 107 and 108, 
local fluctuations eventually collapse, slightly earlier 
in the latter run, confirming earlier speculations by 



extending the results of Hunter et al. (1986) for iso- 
lated gas streams to the fully turbulent regime. In 
fact, given the identical initial conditions and ran- 
dom forcing used in the two runs, the condensation 
that collapses in run 108 also forms in run 107, but is 
unable to actually collapse. 

It is important to emphasize that the collapse in 
the turbulence-induced fluctuations is at small scale, 
contrary to the type of collapse expected for the lin- 
early unstable case. To show this, we have performed 
an additional forced run with J = 1.1 and 7c = 0.9 
(run 109). Contours of the logarithm of the density 
field for runs 109 and 108 are presented in fig. |l^, in 
which the left frame corresponds to run 109 and the 
right frame to run 108. The large void around the 
condensation in run 109 shows that most of the mass 
in the simulation is involved in the collapse, while 
significant amounts of matter are present outside the 
condensation for run 108. Also, note that the "stable" 
run (run 108) takes a rather long time (t = 14) to fi- 
nally develop a collapsing cloud, since a strong enough 
turbulent fluctuation is required, which is statistically 
unlikely. This may be an appropriate mechanism be- 
hind the low efficiency of star formation observed in 
the ISM (see, e.g., [Evans 199i| ). Instead, run 109, 
which is globally unstable, collapses at t = 8.6, a 
time consistent, within a factor of 1.5, with the free- 
fall time iff = Lj/c = 27r/(Jp) = 5.7 in code units. 

Finally, we note that the runs described in this 
section, although providing strong support for the va- 
lidity of the criterion (p7|), cannot be used to deter- 
mine the actual value of 7c, since in order to do this 
it would be necessary to verify that arbitrarily large 
density fluctuations do not collapse, which for large 
7c would require very large Mach numbers. Moreover, 
the dimensionality n of individual shocks may vary, 
and therefore the value of 7c is probably meaningful 
only in a statistical sense. 

4.2. Magnetic field and self-gravity 

We now investigate whether this mechanism is also 
at play in the presence of a uniform magnetic field 
Bq, which we take along the x-direction. The in- 
stability criterion of the uniform-density state is un- 
changed compared to the non-magnetic case (which 
reads J > 1), except for perturbations exactly per- 
pendicular to the field, for which the corresponding 
condition reads > 1 + (Chandrasekhar 1961). 
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Again, note the independence on 7c in these criteria 
in the context of the B-modcl. We have performed a 
forced run similar to run 109 but with _Bo = 2 in code 
units (run 111). This run, with J = 1.1 and 7o = 0.9, 
undergoes contraction along the field lines, but the re- 
sulting slab remains stable in the direction tranverse 
to the field. Note that the critical value Bi of the 
uniform magnetic field for linear stability of the slab 
should be in fact larger than that corresponding to 
the transverse stability of the uniform state. A similar 
run with 7c = 0.3 (run 114) also forms a slab, but the 
turbulence is able to create a subcondensation within 
the slab which collapses gravitationally, thus exhibit- 
ing the same qualitative behavior as that reported 
in the previous subsection for the non-magnetic case. 
This is illustrated in fig. |l^, which shows contour plots 
of the density for both runs. 

The case described above suggests th at the system 
has a supercritical mass to flux ratio (Mouschoviae 
& Spitzer 1976 ), since the magnetic field is not able 
to stop the collapse. This behavior is still observed 
with Bo = 2.5, but at _Bo = 3 the system again 
does not collapse in spite of the production of fluc- 
tuations exceeding p = 55, suggesting that the sys- 
tem has become subcritical. We performed new runs 
with the B-model, no forcing and — '? = 0.04, tak- 
ing as initial conditions zero velocity dispersion and 
the mass concentrated in a wide Gaussian proflle at 
the center of the box with 7e = 0.9, a choice of pa- 
rameter that permits to avoid violent initial collapse 
transients. For both Bo = 'i and Bq = 4 the system 
reaches an equilibrium. Starting now new simulations 
from these equilibrium states with 7c = 0.3 leads, for 
Bo = 4, to a new equilibrium with a higher density 
peak, while for Bo = 3 the system collapses, indi- 
cating that for that value of Bo, thermal pressure is 
still needed to flght gravity. As a result we can argue 
that an arbitrary high level of turbulence would not 
be able to produce gravitational collapse for Bo = 4, 
whereas if a strong enough fluctuation could gather 
enough mass, collapse could occur with Bo = 3. 

Thus, at a fixed value of J, there are two critical 
values of the uniform field Bo', one (i?i) to provide lin- 
ear stability of the slab state, and the other {B2 > Bi) 
corresponding to the transition from supercritical to 
subcritical. According to our results, turbulence is 
able to shift Bi to higher values as turbulent fluc- 
tuations lose thermal support (at small enough 7c). 
However, it does not seem to be able to modify i?2- 
We speculate that Bi should tend to B2 as thermal 



support is weakened (7c 0). Therefore, the mass- 
to-flux ratio criterion appears to be very robust to the 
presence of strong nonlinearities (at the large scale), 
as it is based solely on global conservation conditions. 
However, we note that the presence of small-scale 
MHD turbulence inside the condensations themselves, 
which we cannot follow here due to the limited reso- 
lution, could in principle provide, through reconnec- 
tion, a dissipation of the magnetic fleld allowing for 
an increase of the mass-to-flux ratio. Finally, we re- 
mark that, although turbulence appears incapable of 
inducing collapse of supercritical condensations, it is 
a necessary agent in their formation, since the slab 
state is linearly stable in the absence of turbulence 
due to the contribution from the thermal pressure. 

5. Behavior of a thermodynamic model of the 
ISM at the kpc scale 

In this section we test the applicability of our con- 
clusions on the influence of 7c on cloud formation and 
collapse scenarii in the context of the full model used 
in Paper II (cf. § |^) including also rotation and shear, 
with values adequate for the ISM at the 1 kpc scale. 
We have performed two simulations (runs 102 and 
113) with the flducial values described in Table 1 of 
Paper II, although without star formation nor ran- 
dom forcing (Fg = 0,Fc = 0). These two runs differ 
only in the value of a. We do not include forcing 
since in this case it should be at small scales due to 
stellar activity, which in turn should be the outcome 
of gravitational collapse. 

Neither run 102 nor run 113 collapse gravitation- 
ally, at least over 50 code time units. (Recall that for 
this model, a code time unit corresponds to 1.3 x lO'' 
yr and a density unit equals 1 cm^"^.) Since no forcing 
other than the large scale shear is present in runs 102 
and 113, the turbulence decays after a few turnover 
times and we only observe the nonlinear development 
of the oscillations suggested by Elmegreen (1991) and 
described in Paper II . As in Paper II , the clouds form 
preferentially in regions of zero shear (since shear ro- 
tates a perturbation contracting along the field into 
the perpendicular direction) ; other runs without mag- 
netic field (not shown) form condensations of even 
lower amplitude (due to the absence of magnetic brak- 
ing) and in regions of minimum shear, i.e. a shear 
which minimizes the epicyclic frequency. The absence 
of runaway collapse in these runs is interesting since, 
according to the linear criterion of Paper II, they 
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should be unstable. This confirms the speculation 
by Elmegreen (1991) and the observation of Paper II 
that the development of the instability is not expo- 
nential, but oscillatory, although in the present paper 
the density fluctuations have reached much larger am- 
plitudes, and their ability to rebound is remarkable. 
Ultimately, the possibility that they will collapse on 
very large time scales cannot be ruled out, but is not 
necessarily relevant in the context of the ISM. 

In spite of the absence of gravitational collapse, 
the simulations support the scenario of the previ- 
ous sections: the hard case (a = 0) develops only 
one condensation of fluctuating amplitude, reaching 
Pmax = 20 at t = 14 in code units. Instead, the 
a = 0.5 case develops a larger number of condensa- 
tions, one reaching pmax = 86 at i = 7.6, which how- 
ever rebounds and does not collapse. Later, this same 
condensation reaches Pmax — 95 a.t t — 20.6. The evo- 
lution of the maximum density in the integration box 
is shown for both runs in fig. |l2|. The eventual dying 
out of such fluctuations as seen here for late times 
is linked to the choice Fg — 0, as discussed in Paper 
I and II. These condensations move along epicycles, 
and exhibit a wave-like behavior since there are no 
turbulent compressions nor instabilities of any kind 
(§ 3.1; see also Adams et al., 1994; Gehman et al 



The latter observation is a consequence of the well- 
known fact that the presence of rotation restores 
a "hard" behavior, analogous to a large value of 



7e . I However, as discussed in ? 3.3.1 . the relative 



importance of the Coriolis force compared to the 
magnetic and nonlinear advection terms decreases at 
small scales, so that this scenario suggests stability of 
the large scales with possible instability of the small 
scales. Moreover, the turbulence-induced collapse of 
the small scales shown in § ^ should also contribute 
to small-scale collapse, although in particular in runs 
102 and 113 this effect cannot be seen due to the ab- 
sence of fully developed turbulence and to the low 
resolution. 

6. Conclusions 

In the present paper, we have investigated some 
properties of highly compressible turbulence, which 
is relevant in flows with very efficient cooling pro- 
cesses, such as the ISM or cosmological flows after 
the epoch of recombination. In these cases, the flow 
may be subsonic with respect to the isothermal sound 



speed, and yet have an effectively supersonic behav- 
ior. We first considered the purely hydrodynami- 
cal regime, addressing the type of density structures 
formed by turbulence. We showed that the density 
jump in barotropic flows in the limit of small poly- 
tropic index 7e approaches , where M is the Mach 
number ahead of the shock, rendering the flow much 
more compressible than high Mach number regimes 
with 7e 1. When 7e is small, the density structures 
are long-lived and object-like, i.e. carrying the mass 
together with the density fluctuations while, for larger 
values of 7e, the behavior is more wave- like, and the 
density field consists of short-lived patches separated 
by shocks. 

In the presence of purely compressible forcing, vor- 
ticity is observed to decay in time, even in 3D, where 
the dominant effect of the vortex stretching term is 
actually to act as a sink of potential vorticity at the 
location of shocks. Thus vorticity production relies 
heavily on other sources, such as the Coriolis force, 
large-scale shear originating from differential rotation, 
baroclinic vorticity production in the presence of local 
heating, and the magnetic field. For the Galactic disk, 
the Coriolis force is most important at large scales, 
while the magnetic field is most effective at the small 
ones (molecular clouds). The decay of vorticity we 
observe in the absence of those additional agents may 
be at the origin of the small amounts of rotation found 
in observational surveys of molecular clouds and their 
cores (e.g., Arquilla & Goldsmith 1986; Goodman et 



al. 1993) 



With self-gravity, we have given a simple crite- 
rion for turbulent compression through shocks to pro- 
duce self-gravitating structures, based on the thresh- 
old value of 7c below which the Jeans mass of a turbu- 
lent density fluctuation decreases upon the compres- 
sion. In particular, if the dimensionality of the com- 
pression is between 1 and 2, then the required value of 
7e is between and 1 , consistent w ith estimated val- 
ues in the ISM (e.g., Myers 1978 ). Using numerical 



simulations, we have verified that in the presence of 
strongly compressible turbulence, local turbulent fluc- 
tuations can collapse even if the medium is globally 
stable in the linear sense, although generally requir- 
ing long times. Flows with small 7o develop collapsing 
structures at smaller scales than flows with larger 7e 
but with the same effective Jeans length. Flows with 
self-gravity and magnetic fields exhibit an analogous 
behavior to the purely hydrodynamic case for super- 
critical mass-to-flux ratios. For subcritical situations. 
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large-scale turbulence is incapable of inducing small- 
scale collapse, although it is a necessary agent for the 
formation of magnetically supported clouds. 

These results have implications for present-day vs. 
primordial star formation, inasmuch as 7e is directly 
related to the medium's metallicity. For the low 
metallicity primordial gas, we expect that turbulence 
will not be sufficient to produce small-scale collapsing 
structures due to the large values of 7c, and therefore 
collapse can only occur through large-scale gravita- 
tional instability, leading to the formation of globular 
clusters. Conversely, collapsing turbulence-induced, 
small-scale density fluctuations may lead to formation 
of open clusters under present-day metallicity condi- 
tions. This problem will be investigated in detail in a 
forthcoming paper ( Franco et al. 1996 ). 

Fully thermodynamic models incorporating mag- 
netic fields, rotation and shear with parameters ap- 
prop riate to t he modeling of the ISM at the kpc scale 
(c.f. Paper II ) but without stellar heating (Fg = 0) 
develop highly nonlinear wave-like structures which 
however do not collapse gravitationally. Their be- 
havior is analogous to the hydrodynamic case with a 
larger 7g and initial turbulent transients are not capa- 
ble of forming small-scale self-gravitating structures. 
Thus, in order to obtain collapse of the large scales an 
instability is required. On the other hand, small-scale 
collapse in the clouds should also occur, since rotation 
is less efficient at those scales and we have shown that 
turbulence also induces it. This is not observable in 
the simulations presented in this paper due to their 
low resolution, but may be seen at higher resolutions. 

A final remark is that in this paper we have in- 
vestigated the effects of large-scale turbulence on 
structure formation in highly compressible flows, but 
have neglected small-scale effects, such as those men- 
tioned above, as well as turbulent pressure support 
and reconnection of small-scale magnetic fluctuations, 
which could possibly modify the mass-to-flux ratio 
and allow for an increased compressibility of the flow. 
The simultaneous investigation of both aspects of 
turbulence requires very high resolution simulations, 
which will be presented elsewhere. 
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Fig. 1. — Contours of the density field for run 63 
(7c — 1.5, left) and run 65 (70 — 0.3, right) at a 
nondimensional time oi t = 7.6. In both cases there 
are 10 contour levels with a constant spacing Ap = 
0.3. Note the "patchy" structure of the density field, 
with higher peaks at the intersection of shocks in (b) . 



Fig. 2. — Contours of the density field of run 54 
{a — —0.5, left) and run 53 (a = 0.5, right) at 
t = 6.0. The contours represent density values of 
p = 1,2,4,8,16,32. Both runs use the ppm equa- 
tions. Note that clouds are better defined than in fig. 
|l|, and have higher peak densities in (b). 
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Fig. 3. — One-dimensional sketches of the density p, 
the passive scalar s and the velocity u in the case of 
a wave (top) and an object (bottom). In the case of 
the wave, the passive scalar, which represents advec- 
tion of material, drifts away from the density, while 
in the case of the object, s remains within the density 
fluctuation. 
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Fig. 4. — Evolution of the ratio of solenoidal to total 

kinetic energy per nnit mass Cg/ek for runs 60 (solid 
line), run 67 {dotted line), run 68 {short-dashed line) 
and run 69 {long-dashed line). Runs 60 and 67 are 
in 3D and use the ppm equations with a = 0.5 and 
a = —0.5 respectively. Run 68 is similar to run 60 
but fully thermodynamic. Run 69 is similar to run 60 
but includes the Coriolis force with = 0.4. After 
similar initial transients, the ratio Cg/ek decays in all 
cases except in the run including the Coriolis force. 




Fig. 5. — Gray-scale images of the density field for 
run 28 {left) and run 35 {right). Both runs are fully 
thermodynamic, but run 28 has local thermal forcing 
only, while run 35 has wind-type compressible forcing 
on the velocity only. The shells evolve smoothly in 
run 35, but distort and disrupt in run 28 due to the 
strong vorticity production. 
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Fig. 6. — Evolution of the ratio of solenoidal to to- 
tal kinetic energy per unit mass Cg/ek for runs 86 
{solid line), run 84 {dotted line), run 87 {short-dashed 
line) and run 88 {long-dashed line). These runs re- 
spectively have values of the initial uniform field of 
Bo = 0.05, 0.3, 1 and 3. The final value of e^/e-k is seen 
to increase with Bq until saturation at the equiparti- 
tion level is reached for Bo > 1. 
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Fig. 7. — Evolution of the fluctuating magnetic en- 
ergy for the same runs as in fig. ^, using the same hne 
coding. A similar trend of increasing strength with 
Bo is observed, with saturation at > 1. Note also 
the increasing frequency of the oscillations, roughly 
proportional to B^- 



Fig. 9. — Evolution of the maximum density for runs 
106 {solid line), 107 (dotted line) and 108 (dashed 
line), which respectively have 70 = 0.9,0.3 and 0.1, 
all with J = 0.9, implying linear stability. Turbulent 
fluctuations of flnite strength are seen to be able to 
collapse gravitationally for runs 107 and 108 in spite 
of the global linear stability. 
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Fig. 8. — Spectra of the x- and y-components of 
the magnetic field of run 87 at i = 14.4. A well- 
defined power-law inertial range is observed, indica- 
tive of fully developed turbulence. 



Fig. 10. — Contour plots of the density fields of runs 
108 (left) and 109 (right), which respectively have J = 
0.9 and 1.1. Run 109 is thus linearly unstable, and 
the collapse is seen to involve most of the mass in 
the simulation, while the turbulence-induced collapse 
in run 108 is local. The contours denote the values 
/9 = 0.5, 1,2, 4, 8, 16, 32 and 64. 
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Fig. 11. — Contour plots of the density fields of runs 
111 {left) and 114 (right), both permeated by a uni- 
form magnetic field of strength Bo = 0.2 in the x- 
direction and with J = 1.1. Run 111 has 7e = 0.9, 
and for run 114, 7c — 0.3. Being linearly unstable, 
both runs contract to a slab. Turbulent fluctuations, 
however, manage to produce a subcondensation in run 
114 which finally collapses gravitationally. The con- 
tour spacing is identical to that of fig. nol 




Fig. 12. — Evolution of the maximum density for runs 
102 {solid line) and run 113 {dotted line). Both runs 
include rotation, shear and magnetic field, in the con- 
text of the full model of Paper II, but without star 
formation. Although not collapsing gravitationally, 
both runs exhibit strong density fluctuations. 
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Table 1. Run parameters 



Run # 


Type" 


Forcing'' 


(e,/et.)n'^ 


Shear'' 


J"" 


bJ 

u 




7/ 

/ e 


a' 


28 


T2D 


W 





N 








0.0 


N/A 


0.5 


35 


T2D 


T 





N 








0.0 


N/A 


0.5 


53 


ppni2D 


R 




N 











N/A 


0.5 


54 


ppni2D 


R 


I 


N 











N/A 


-0.5 


57 


ppm2D 


R 


I 


N 








0.4 


N/A 


0.5 


59 


ppni2D 


R 


I 


Y 








0.0 


N/A 


0.5 


60 


ppmSD 


R 


I 


N 











N/A 


0.5 


63 


B2D 


R 


I 


N 











1.5 


N/A 


65 


B2D 


R 


^ 


N 











.3 


N/A 


67 


ppmSD 


R 




N 











N/A 


-0.5 


68 


T3D 


R 




N 











N/A 


0.5 


69 


ppmSD 


R 




N 








0.4 


N/A 


0.5 


84 


ppni2D 


R 




N 





0.3 





N/A 


0.5 


85 


ppm2D 


R 


1 


N 








0.05 


N/A 


0.5 


86 


ppm2D 


R 




N 





0.05 





N/A 


0.5 


87 


ppni2D 


R 




N 





1. 





N/A 


0.5 


88 


ppm2D 


R 




N 





3. 





N/A 


0.5 


102 


T2D 


N/A 




Y 


.51 


.3 


0.4 


N/A 


0.5 


106 


B2D 


R 




N 


0.9 


0.0 





0.9 


N/A 


107 


B2D 


R 




N 


0.9 


0.0 





0.3 


N/A 


108 


B2D 


R 




N 


0.9 


0.0 





0.1 


N/A 


109 


B2D 


R 




N 


1.1 


0.0 





0.9 


N/A 


111 


B2D 


R 




N 


1.1 


2.0 





0.9 


N/A 


113 


T2D 


N/A 




Y 


.51 


.3 


0.4 


N/A 


0. 


114 


B2D 


R 




N 


1.1 


2.0 





0.3 


N/A 
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Table 1 

Main characteristics of the runs described in this 
paper, together with the run number. AU runs except 
run 35 exclude stellar heating (i.e. with Fg = in (j^). 
Computations are done on a uniform grid of either 
128^ or 64^ points, using a pseudo-spectral method 
with hyperviscosity. For more detail on the imple- 
mentation of the physical terms and in the nondi- 
mensionalization of the equations, see text and Paper 
I and IF 

° Types of computations: T2D (resp. T3D) refers 
to the full thermodynamical equations (|l|) - (^) writ- 
ten in §2, integrated in two space dimensions (resp. 
three); ppm2D (resp. ppmSD) refer to a simplifica- 
tion of the full equations in which a piecewise poly- 
tropic model (see (p^)) is used in two (resp. three) 
dimensions; and finally B2D (resp. BSD) correspond 
to the further simplification of a barotropic law in two 
(resp. three) dimensions. 

* The type of forcing used in the different models can 
be either: (R), a random forcing with a steep (^ fc~'*) 
spectrum; (T), a thermal energy forcing with Fg ^ 0; 
or finally (W), a wind with Fc = —Vi/j with ^ ~ Fs. 

(es/efc)o is the ratio at t = of the rotational to 
total kinetic energy. 

Y (resp. N) indicates the presence (resp. absence) 
of shear in the computation. 

J is the gravitational parameter and measures the 
number of Jeans masses in the flow; when J = 0, no 
self-gravity is included. 

^ Bo is the mean magnetic field; Bo = 1 corresponds 
to an Alfven velocity equal to the rms turbulent ve- 
locity (see Paper II). In this column, Bo = indicates 
a non-magnetic run. 

s ^l measures the intensity of the Coriolis force, in- 
dependently of the externally imposed shear. 

^ 7e is the effective polytropic exponent used in the 
purely barotropic (B) models (see text). 

* a is the exponent of the density used in a model 
of diffuse heating (see Paper II) . The relationship be- 
tween a and the set of effective polytropic ex- 
ponents appearing in the ppm models is given in eq. 
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